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ABSTRACT 

We describe the data reduction algorithms for the Multiband Imaging Photometer for Spitzer 
(MIPS) instrument. These algorithms were based on extensive preflight testing and modeling of 
the Si:As (24 /jm) and Ge:Ga (70 and 160 /im) arrays in MIPS and have been refined based on initial 
flight data. The behaviors we describe are typical of state-of-the-art infrared focal planes operated in 
the low backgrounds of space. The Ge arrays are bulk photoconductors and therefore show a variety 
of artifacts that must be removed to calibrate the data. The Si array, while better behaved than the 
Ge arrays, does show a handful of artifacts that also must be removed to calibrate the data. The 
data reduction to remove these effects is divided into three parts. The first part converts the non- 
destructively read data ramps into slopes while removing artifacts with time constants of the order of 
the exposure time. The second part calibrates the slope measurements while removing artifacts with 
time constants longer than the exposure time. The third part uses the redundancy inherit in the MIPS 
observing modes to improve the artifact removal iteratively. For each of these steps, we illustrate the 
relevant laboratory experiments or theoretical arguments along with the mathematical approaches 
taken to calibrate the data. Finally, we describe how these preflight algorithms have performed on 
actual flight data. 

Subject headings: instrumentation: detectors 



1. INTRODUCTION 

Most of our knowledge of the Universe at far infrared 
wavelengths has been obtained with photoconductive de- 
tectors, particularly as used in the Infrared Astronomy 
Satellite (IRAS) and the Infrared Space Observatory 
(ISO). These detectors have been selected because they 
provide excellent performance at relatively elevated op- 
erating temperatures (compared with those needed to 
suppress thermal noise in bolometers) . Similar consider- 
ations led to development of high performance photocon- 
ductor arrays for the Multiband Imaging Photometer for 
Spitzer (MIPS), namely a Ge:Ga array and a stressed 
Ge:Ga array operating at 70 and 160 /im respectively 
which were built at the University of Arizona. To provide 
complementary measurements at 24 /im, the instrument 
also includes a Si:As Blocked Impurity Band (BIB) array, 
built at Boeing North America (BNA) under contract to 
the Infrared Spectrograph (IRS) team. 

In the BIB or Impurity Band Conduction (IBC) archi- 
tecture, the high impedance required to minimize John- 
son noise is provided by a thin, high-purity layer of sil- 
icon. The infrared absorption occurs in a second layer, 
which can be relatively strongly doped. Due to the sep- 
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aration of these two functions, the detector layers can 
be optimized separately. Thus, these devices can be de- 
signed and built to have fast response, high resistance to 
cosmic ray irradiation induced responsivity shifts, high 
quantum efficiency, and good photometric behavior. Be- 
cause the processing necessary for high performance sili- 
con IBC devices has only relatively recently become pos- 
sible there is relatively little experience with them in 
space astronomy missions. An early-generation detec- 
tor array was used in the Short Wave length Spectrome- 
ter (SWS) in ISO ijKessler et alJlpjol . Initially, the de- 
vice showed degradation due to damage by large ionizing 
particle exposures when the satellite passed through the 
trapped radiation belts. Once the operating conditions 
were adjusted to minimize these effects, the SWS detec- 
tors showed the expected virtues of this type of device 
even though they d id not achieve their preflight sensi- 
tivity expectations llde Oraauw et all 119961 : iHeras et alJ 
2000l lValentiin fc Thill2n0flL 

At wavelengths longer than 40 /im, photoconductors 
are typically built in germanium because of the avail- 
ability of impurity levels in this material that are much 
more shallow than those in silicon. Achieving the appro- 
priate structure and simultaneously the stringent impu- 
rity control for germanium IBC devices has proven dif- 
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ficult. As a result, all far infrared space astronomy mis- 
sions have used simple bulk photoconductors. MIPS uses 
bulk gallium doped germanium (Ge:Ga) detectors, both 
unstressed and stressed. In such detectors, the same vol- 
ume of material determines both the electrical and photo- 
absorptive properties, making the optimization less flex- 
ible than with Si IBC detectors. Consequently, their be- 
havior retains some undesirable properties that can be 
circumvented through the more complex architecture of 
IBC devices. Nonetheless, generally satisfactory perfor- 
mance is possible and has been achieved in past space 
astronomy missions. 

Th e 60 and 100 /jm bands in IRAS ijNeugebauer et all 
1984) utilized 15 Ge:Ga photoconductors each. The de- 
tectors were read out with transimpedance amplifiers 
that used junction field effect transistor (JFET) first 
stages mounted in a way that isolated them thermally. 
This allowed these transistors to be heated, resulting 
in low noise an d stable operation ijRieke et alJ 119811: 
ILow et aX1ll984|) . Detector calibration was maintained 
by flashing reverse bolometer stimulators mounted in the 
center of the telescope secondary mirror, and cosmic ray 
effects were erased by boostin g the detector bias to break- 
down (Be ichman et al.l|l988|) . The intrinsic performance 
of the detectors was limited by the Johnson noise of the 
transimpedance amplifier (TIA) feedback resistors and 
by other noise sources associated with the readout. The 
in-flight performance was similar to expectations from 
pre-flight calibrations. 

The ISOPHOT instrument l|Lemke et al.ll996ft carried 
a 3 x 3 array of unstressed Ge:Ga detectors operating 
from 50 to 105 /im and a 2 x 2 array of stressed de- 
vices operating from 120 to 200 /.im. The readout was by 
a capacitive transimpedance metal oxide semiconductor 
field-effect transistor (MOSFET) amplifier whose pro- 
cessing had been adj usted to improve it s performance 
at low temperatures l|Dierckx et al.lfl"992^ . Calibration 
was assisted with a stimulator built into the instrument, 
which could be viewed by adjusting the position of a 
scan mirror. In practice, the unstressed focal plane never 
achieved the performance level anticipated from labora- 
tory measurements of its noise equivalent power (NEP). 
The performance of the stressed devices was substan- 
tially better, due in part to the relatively large fast re- 
sponse component of these devices (compared with the 
slow component) and their better thermal isolation from 
the readout amplifiers. 

The L ong Wavelength Spect rometer (LWS) instrument 
on ISO ijSwinvard et al.lll996j) used a single Ge:Be detec- 
tor, five Ge:Ga detectors, and four stressed Ge:Ga detec- 
tors. The readouts were based on JFETs, mounted with 
thermal isolation and heated to a temperature where 
they operated with good stability and low noise. The 
readout circuit was an integrating source follower and 
NEPs of ~ 1 x 10~ 18 W/Hz 1 / 2 were measured in the 
laboratory l|Church et alJl993T) . Calibration was assisted 
with built in stimulators that were flashed between spec- 
tral scans. On orbit, it was found that frequent small 
glitches, probably associated with cosmic ray hits, lim- 
ited the maximum integration times to shorter values 
than had been a nticipated and also required a lower 
opera ting voltage ijBuredorf et al.lll998t iSwinvard et al.l 
2000). With these mitigations, the NEPs were found to 



be ^ 4 higher in orbit tha n expected from ground test 
data | ISwinvard et aDISOQCh . 

In MIPS, the Ge:Ga detectors are carefully isolated 
thermally from their readouts and operated at suffi- 
ciently cold temperatures that their dark currents are 
low and stable. The MOSFET-based readouts use a spe- 
cialized foundry process that provides them with good 
DC stability even at the low operating temperature of 
~1.5K. This feature, combined with the capacitive TIA 
(CTIA) circuit, maintains the detector bias accurately. 
A scan mirror (based on a design provided by T. dc 
Graauw) modulates the signals on a pixel so measure- 
ments can be obtai ned from the relatively well-behaved 
ijHaegel et al.ll200ll) fast component of the detector re- 
sponse. Responsivity variations are tracked with the aid 
of frequent stimulator flashes. Finally, the instrument 
operations force observers to combine many short obser- 
vations of a source into a single measurement. The high 
level of redundancy in the data helps identify outlier sig- 
nals and also improves the calibration by simple aver- 
aging over variations. The efficacy of this operational 
approach is confirmed by the on-orbit results. Details 
o n the design and construction of MIPS c an be found 
inlHeim et alJ jl998|) : ISchnurr et all l)1998ft : lYoung et alJ 
(1998) . The inflight performance of MIPS is described 
bv lRieke et all (|2004tL 

This paper describes the approaches for reduction and 
calibration of the MIPS data. Section [3 details the chal- 
lenges of using Si and Ge detectors in a space astronomy 
mission. Section|3gives a summary of the design and op- 
erational features of MIPS that address these challenges. 
Section 0] gives an overview of the three stages of MIPS 
data reduction. These stages are discussed in more detail 
in the following three sections. Sectional details the pro- 
cessing steps to turn the integration ramps into measured 
slopes. Section discusses the corrections to transform 
the slopes into calibrated fluxes. Section gives a brief 
overview of the use of the inherent redundancy in the 
observations to further improve the reduction. Section|H] 
gives the results of initial testing of these reduction tech- 
niques with flight data. Finally, Section El provides a 
summary. 

2. THE CHALLENGE 
2.1. Germanium Arrays (70 & 160 /im) 

At high backgrounds, such as might be encountered 
in an airborne instrument, far infrared photoconductors 
behave relatively well, with rapid adjustment of the de- 
tector resistance appropriate to a change in illumination 
level. As the background is decreased, the adjustment to 
equilibrium levels occurs in a multistep process with mul- 
tiple associated time constants as discussed below. Thus, 
the detectors can be used in a straightforward manner at 
high backgrounds but precautions must be taken at low 
ones to track the cali bration. For a more detailed discus- 
sion see lRiekd (j2002lL 

The fast response component in these detectors results 
from the current conducted within the detector volume 
associated with the drift of charge carriers freed by ab- 
sorption of photons. The speed of this component is 
controlled by the propagation of a zone boundary with 
drift velocity v^, so that the time constant is given by 
the free carrier lifetime divided by the photoconductive 
gain. This time is very fast (microseconds or shorter) in 
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comparison to normal detection standards. However, as 
charge moves within the detector, the electrical equilib- 
rium must be maintained. For example, charge carriers 
generated by photoionization are removed from the de- 
tector when they drift to a contact. They are replaced by 
injection of new charge carriers from the opposite con- 
tact, but the necessity for new charge can only be con- 
veyed across the detector at a characteristic time pro- 
portional to the "dielectric relaxation time" , basically 
its capacitive or RC time constant: 



Here, kq is the dielectric constant of the material and 
/i is the mobility for the charge carrier of interest, e is 
the permittivity of free space, Uq is the density of free 
carriers, and q is the charge of the electron. 

The slow response components arise from this phe- 
nomenon. The form of this time constant makes explicit 
the dependence on illumination level through the density 
of free charge carriers, n a . In fully illuminated detec- 
tors (for example, the integrating cavities used for the 
160 [im array) and at the low backgrounds appropriate 
for space-borne operation, can be tens of seconds. In 
transverse contact detectors, such as those used for the 
MIPS 70 /jm array, the part of the detector volume near 
the injecting contact may be poorly illuminated and have 
large resistance. The detector therefore adjusts to a new 
equilibrium only at the large dielectric time constant of 
this layer, which can be hundreds of seconds at low back- 
grounds. The initial shift of charge in the detector can 
set up a space charge that reduces the field in the bulk 
of the device, leading to a reduction of responsivity fol- 
lowing the initial fast response. From its appearance on 
a plot of response versus time, this behavior is described 
as "hook" response. As the field is restored at a char- 
acteristic rate of r<j, th e response grows slo wly to a new 
equilibrium value. See lHaegel et alJ (|2001^ for detailed 
modeling of these effects. 

In space applications, ionizing particles such as cosmic 
rays also affect the calibration of these detectors. The 
electrons freed by a cosmic ray hit can be captured by 
ionized minority impurities, reducing the effective com- 
pensation and increasing the responsivity. The shifts in 
detector characteristics can be removed by warming it to 
a temperature that re-establishes thermal equilibrium, 
and then cooling it back to proper operating conditions. 
Between such anneal cycles, the responsivity needs to 
be tracked to yield calibrated data. All successful uses of 
far infrared photoconductors at low backgrounds have in- 
cluded local relative calibrators of reverse bolometer de- 
sign that allow an accurately repeatable amount of light 
to be put on the detector. These stimulators allow fre- 
quent measurement of the relative detector responsivity. 
In general, this strategy is most successful when the con- 
ditions of measurement are changed the least to carry out 
the relative calibration. The MIPS instrument includes 
such calibrators, which are flashed approximately every 
two minutes. Based upon data obtained at a proton ac- 
celerator and in space, the average increase in response 
over a two minute period in the space environment can 
be 0.5% to 1%, so the calibration interval allows tracking 
the response accurately. 



2.2. Silicon Array (24 urn) 

Although the detectors in the silicon array are expected 
to perform well photometrically, the array as a system 
shows a number of effects that must be removed to ob- 
tain calibrated data. The array is operated well below 
the freezeout temperature for the dopants in the silicon 
readout (the readout circuit uses a different foundry pro- 
cess from that developed for the Ge detectors). There- 
fore, the array must be operated in a continuous read 
mode to avoid setting up drifts in the outputs that would 
degrade the read noise. The flight electronics and soft- 
ware are designed to maintain a steady read rate of once 
per half MIPS second (see i]3.5(l . When the array is first 
turned on, the transient effects of the readout cause a 
slow drift in the outputs. Much of this effect can be 
removed by annealing the array, which is the standard 
procedure for starting the MIPS 24 (im operations. 

The array shows an effect termed "droop." The output 
of the device is proportional to the signal it has collected, 
plus a second term that is proportional to the average sig- 
nal over the entire array. In addition, the 24 /jm array 
has a number of smaller effects (e.g., rowdroop, electronic 
nonlinearities, etc.) which are described later in this pa- 
per. 

3. DESIGN AND OPERATION OF MIPS 

The design and operation of MIPS is summarized here, 
paying special attention to those areas which answer the 
challenges outlined above and, therefore, produce data 
that can be reduced successfully. 

3.1. Instrument Overview 

MIPS has three instrument sections, one for 24 /im 
imaging, one for 70 /im imaging and low resolution spec- 
troscopy, and one for 160 /im imaging. Light is directed 
into the three sections off a single axis scan mirror. 

The 24 /im section uses a 128 x 128 pixel Si:As IBC 
array and operates in a fixed broad spectral band ex- 
tending from 21 to about 27 /jm (the long wavelength 
cutoff is determined by the photo-absorptive cutoff of 
the detector array). After light enters this arm of the 
instrument from a pickoff mirror, it is brought to a pupil 
on a facet of the scan mirror. It is reflected off this mirror 
into imaging optics that relay the telescope focal plane to 
the detector array at a scale of 2'.'5 per pixel correspond- 
ing to a X/2.2D sampling of the point spread function, 
where D is the telescope aperture. The surface area of a 
24 /im pixel is 75 x 75 /im 2 . The field of view provided 
by this array is 5'3. A reverse bolometer stimulator in 
this optical train allows relative calibration signals to be 
projected onto the array. The scan mirror allows images 
to be dithered on the array without the overheads as- 
sociated with moving and stabilizing the spacecraft. It 
also enables an efficient mode of mapping (scan mapping) 
in which the spacecraft is scanned slowly across the sky 
and the scan mirror is driven in a sawtooth waveform 
that counters the spacecraft motion, freezing the images 
on the detector array during integrations. 

The 70 /im section uses a 32 x 32 pixel Ge:Ga array 
sensitive from 53 to 107 /^m. A cable failure external to 
the instrument has disabled half of the array and the fol- 
lowing description reflects this situation. The light from 
the telescope is reflected into the instrument off a second 
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pickoff mirror. It is brought to a pupil at a second facet 
of the scan mirror and from there passes through optics 
that bring it to the detector array. For this arm of the 
instrument, there are actually three optical trains that 
can relay the light to the array; the scan mirror is used 
to select the path to be used for an observation. One 
train provides imaging over a field 2.'7 x 5.'3, with a pixel 
scale of 9'.'8 corresponding to a A/1.81? sampling of the 
point spread function. The physical size of a 70 /jm pixel 
is 0.75 x 0.75 mm 2 and 2 mm long in the direction of 
the optical axis. This train provides imaging over a fixed 
photometric band from 55 to 86 /im. The scan mirror 
feeds this mode when it is in position to feed the other 
two arrays, so imaging can be done on all three arrays 
simultaneously. A second train also provides imaging in 
the same band, but with the focal plane magnified by a 
factor of two to 4'.'9 per pixel. This mode is provided 
for imaging compact sources where the maximum possi- 
ble angular resolution is desired: the pixel scale corre- 
sponds to A/3.5D at the center wavelength of the filter 
band. The third train brings the light into a spectrom- 
eter, with spectral resolution of R = A/AA ^25—15 
from 53 — 107 /jm. In this Spectral Energy Distribution 
(SED) instrument mode, light is directed to a reflective 
"slit" and then to a concave reflective diffraction grating 
that disperses the light and images the spectrum onto a 
portion of the 70 /im array. The slit is 16 pixels long 
and 2 pixels wide, corresponding to 2.'7 x 0.'32 on the 
sky. The dispersion is 1.73 (im pixel -1 . Reverse bolome- 
ter stimulators are provided for calibration, and the scan 
mirror provides the dithered and scan mapping modes of 
operation at 70 /jm as have been described for the 24 /im 
array. 

The 160 /im section shares the pickoff mirror and scan 
mirror facet with the 24 /im band. After the light has 
been reflected off the scan mirror, the telescope focal 
plane is reimaged and divided, with part going to the 
Si:As array and part going to a stressed Ge:Ga array, op- 
erating in a fixed filter band from 140 to 180 /im. This 
array has 2 x 20 pixels, arranged to provide an imaging 
field 5.'3 long in the direction orthogonal to the scan mir- 
ror motion with the two rows of detectors spaced such 
that there is a gap one pixel wide between the two rows. 
This pixel size provides X/2.2D sampling of the point 
spread function. The physical size of a 160 /im pixel is 
0.81 x 0.81 x 0.81 mm 3 . Reverse bolometer stimulators 
are included in the optical train, and the scan mirror pro- 
vides modes similar to those with the other two arrays. 

3.2. Stimulators 

A key aspect of the calibration of the MIPS Ge ar- 
rays is the frequent use of stimulators ijBeeman fc Hailed 
2002) to track responsivity variations. The emitters in 
these devices are sapphire plates blackened with a thin 
deposition of bismuth, which also acts as an electrical 
resistor. The emitters are suspended in a metal ring by 
nylon supports and their electrical leads. When a con- 
trolled current is run through the device, the sapphire 
plate is rapidly heated by ohmic losses in the metallized 
layer. The thermal emission is used to track changes in 
detector response in a relative manner; hence these de- 
vices are described as stimulators rather than calibrators. 
Because of the large responsivity of the detector arrays, 
it is necessary to operate these stimulators highly inef- 



ficiently to ensure accurate control without blinding the 
detectors. They are mounted inside cavities that are in- 
tentionally designed to be inefficient (e.g., black walls, 
small exit holes). This allows the stimulators to be run 
at high enough voltage to be stable and emit at a rea- 
sonable effective temperature. 

The constant-amplitude stimulator flashes provide a 
means of tracking the responsivity drift inherent in the 
Ge detectors. Figure Q illustrates the importance of 
tracking the responsivity variations of Ge detectors with 
as fine a time resolution as feasible. The repeatability 
of the stimulator measurements is a function of both the 
background seen by the detector element as well as the 
amplitude of the stimulator (stim) signal above the back- 
ground. The repeatability of a measurement of the stim 
signal improves with decreasing background and increas- 
ing stim amplitude. For both the 70 and 160 /im ar- 
rays, in ground testing stim amplitudes of greater than 
~7500 DN/s above the background yielded a repeatabil- 
ity of better than ~1% on most backgrounds. Setting 
the stim amplitudes at this level provides a balance be- 
tween repeatability of the stims and the range of back- 
grounds accessible to observation without saturation. At 
this level, well over 95% of the sky should be observable 
without saturating stim flash measurements at both 70 
and 160 /im. 

Additional complications at 160 /im include a strong 
illumination gradient in the stim flash illumination pat- 
tern from one end of the array to the other as well as a 
large increase in the responsivity of the array with ex- 
posure to cosmic rays. It is not possible to set the stim 
amplitude at the optimum 7500 DN/s across the whole 
array due to a factor of four gradient in the stim ampli- 
tude across the array. The on-orbit stim amplitude was 
set to provide an optimal amplitude over the majority of 
the 160 /im pixels. The degradation in stim repeatabil- 
ity on the low illumination region can be mitigated by 
an observing strategy that dithers the image such that 
the same region of the sky spends equal amounts of time 
on both regions of the detector. 

3.3. Anneals 

Both of the Ge arrays show calibration shifts with even 
small exposure to ionizing radiation. The effects of ion- 
izing particles were tested using characterization arrays 
(see m at the University of California, Davis accel- 
erator. The proton beam was attenuated to reduce the 
particle impact rate to a level similar to that expected 
on orbit. The energy of the particles was such that each 
impact was strongly ionizing, depositing much more en- 
ergy in the detector volume than is expected from a typ- 
ical cosmic ray. Thus, these tests served as a worst-case 
model of the detector response to cosmic-rays on orbit. 

The detector responsivity slowly increased with time 
under exposure to the proton beam. The rate of respon- 
sivity increase on the 70 /im array was comparable to 
that observed under typical illumination conditions (cf. 
Fig. ^) without the proton beam, suggesting that ac- 
cumulated transient response from the background and 
signals inside the cold test chamber and the photon flux 
at the accelerator contribute similarly to the responsivity 
increase. If the particle impacts at the accelerator really 
represent a worst-case scenario, this suggests that the 
on-orbit responsivity increase of the 70 /im array may be 
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Fig. 1. — The stim flash amplitudes for a single pixel of the 70 fim (left) and 160 fim (right) arrays are plotted from a 2 hour test where 
the stim was flashed approximately every 2 minutes on a constant background. The vertical dotted lines denote the image in which a 
cosmic ray was detected. On the 70 fim plot, the first cosmic ray can be seen to have caused a ~3%responsivity increase. The unit DN/s 
refers to Data Numbers per second. The x-axis gives the number of images (DCE is defined in H3.5I taken where each image represents 
about 1 1 seconds of data. 



dominated by photon flux rather than cosmic ray effects. 
In contrast, the 160 fim array showed a large responsivity 
increase with increasing radiation dose. 

If they are of modest size, such responsivity shifts can 
be determined and removed during calibration through 
use of the stimulator observations. However, when the 
shifts are large, they are also highly unstable and can 
result in substantial excess noise. Three methods were 
tested to remove such effects: rc-thermalization of the 
detectors by heating them (anneals) , exposing the detec- 
tors to a bright photon source, and boosting the detector 
bias above breakdown. Our experiments indicated that 
the latter two methods produced little benefit. Although 
the instrument design permits use of all three techniques, 
we remove radiation damage to the Ge arrays by period- 
ically thermally annealing the detectors. 

3.4. Observing Modes 

There are four MIPS observing modes, all of which 
have been designed to provide a high level of redundancy 
to ensure good quality data (especially for the Ge arrays). 
The Photometry mode is for point and small sources. As 
an example of the redundancy inherit in MIPS observa- 
tions, a visualization of a single Photometry Mode cy- 
cle is shown in Fig. [3 for 70 /im. The Scan Map Mode 
provides efficient, simultaneous mapping at 24, 70, and 
160 /im by using a ramp motion for the scan mirror to 
compensate for continuous telescope motion, effectively 
freezing images of the sky on the arrays. A visualization 
of a small portion of a scan leg is shown in Fig. |21 The 
SED mode provides 53 to 107 /jm spectra with a resolu- 
tion R 25 - 15. The Total Power Mode (TPM) is for 
making absolute measurements of extended emissions. 

3.5. Data Collection 

The pacing of the MIPS data collection is based on a 
"MIPS second." A MIPS second is approximately 1.049 
seconds, and has been selected to synchronize the data 
collection with potential sources of periodic noise, such 
as the computer clock or the oscillators in the power 
supplies. To first order, this design prevents the down- 
conversion of pickup from these potential noise sources 
into the astronomical signals. The data are taken in Data 



Collection Events (DCEs); at the end of a DCE, the ar- 
ray is reset before taking more data. DCEs are currently 
limited to 3, 4, 10, or 30 MIPS seconds for the 24 fim ar- 
ray and 3, 10, or 10 MIPS seconds for the 70 and 160 fim 
arrays. 

During a DCE, each pixel generates a voltage ramp on 
the array output, as the charge from incoming photons is 
accumulated on the input node of its integrating ampli- 
fier. These ramps are the basic data collected by all three 
arrays. The 24 /im array is non-destructively read out 
every 1/2 MIPS second while the 70 and 160 fim arrays 
are non-destructively read out every 1/8 MIPS second. 
All the samples are downlinked for the 70 and 160 fim 
arrays, but this is not possible for the 24 /im array due to 
bandwidth restrictions. The 24 fim array has two data 
modes, SUR and RAW. Most 24 fim data are taken in 
SUR mode in which the ramps are fitted to a line, and 
only the fitted slope and first difference (the difference 
between the first two reads in the ramp) are downlinked. 
The RAW mode downlinks the full 24 fim ramps, but 
this mode is used only for engineering observations. 

4. OVERVIEW OF MIPS DATA PROCESSING 

There are three natural steps in reducing data from 
integrating amplifiers: (1) converting the integration 
ramps to slopes; (2) further time-domain processing of 
the slope images; and (3) processing of dithered images in 
the spatial domain. For detectors that do not have time- 
dependent responsivities, only the first and last steps are 
usually important. This is strongly not the case for the 
MIPS Ge arrays and also mildly not so for the MIPS Si 
array. 

As a result, MIPS processing includes all three steps. 
First, the integration ramps are converted into slopes 
(DN/s) while removing instrumental signatures with 
time constants on the order of the DCE exposure times 
(iJHJ. Second, the slopes are calibrated and instrumental 
signatures with time constants longer than the DCE ex- 
posure times are removed (®. Third, the redundancy 
inherent in the MIPS observing modes allows a second 
pass at removing instrumental signatures (®. The al- 
gorithms used in the first two steps have mostly been 
determined. The main algorithms used by the third step 
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Fig. 2. — The images for a single observation cycle in 70 fim Compact Source Photometry mode. The numeral on each individual image 
gives the object image number in the cycle. The "FLASH" designation corresponds to a stim flash image and the image taken before the 
stim flash gives the background on top of which the stim is flashed. The "COMBfNED" image was made with pixels 1/4 the original size. 
Note that an object in the center of the combined image is significantly better sampled than one near the edges. 




24^.m mosaic 



Fig. 3. — A visualization for the Scan Map mode. The stim flashes are not shown. On the left, the locations of 5 simultaneous images at 
24, 70, and 160 £im are shown on a sky image at the 24 (im resolution. The single field of view for each array is denoted by a bold outline. 
The middle shows the individual images at 24 and 160 fivci; the 70 /im images resemble the 24 (im images, just offset downward. On the 
right, the 24 and 160 fira mosaics created from the 5 individual frames are shown. 



are being optimized with actual data taken on orbit. Por- 
tions of the reduction algorithms descri bed in this paper 
were p resented in a preliminary form by Hcsscl roth et all 
(200(|). 

We made extensive use of laboratory testing and the- 
oretical investigations in choosing and ordering the rele- 



vant steps. Figure El is a graphical representation of the 
specific tasks in each of the three processing steps. 

4.1. Laboratory Testing of Ge Arrays 

Three versions of the 70 and 160 /im arrays were con- 
structed: a flight array, a flight spare array, and a charac- 
terization array. Before integration into the instrument, 
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Fig. 4. — Graphical representation of the flow of the reduction of MIPS data. 



the basic performance of the flight and flight spare ar- 
rays was measured (i.e., read noise, dark current, NEP, 
etc.). The characterization arrays were then installed in 
the two specialized dewars previously used for the flight 
and flight spare array testing. These arrays are used to 
determine the detailed behaviors of the 70 and 160 /im 
detectors. This knowledge was then used to design ob- 
servations with the flight arrays to remove specific Ge 
detector effects. The ability to do extensive testing on 
the characterization arrays has been crucial to the de- 
velopment of the data reduction algorithms for the Ge 
arrays detailed in this paper. 
In addition to testing at the array level, testing at the 



instrument level was carried out using the Low Back- 
ground Test Chamber (LBTC). The LBTC was con- 
structed to allow for testing of the full MIPS instrument 
and, thus, had a number of independently controlled 
stimulators including pinhole stimulators providing point 
sources for testing. The LBTC allowed for the imaging 
performance of the full instrument to be tested as well 
as providing for extensive testing of the 24 /im array. 

Additi onal details of the laboratory testing can be 
found in lYoung et af] l|2003|) . 



4.2. Numerical Modeling of Ge Arrays 
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read § 

Fig. 5. — Example of a 24 /an RAW data ramp for a 30 second 
DCE, showing several steps along the processing pipeline. Aster- 
isks represent the ramp after dark subtraction and show ADC sat- 
uration is reached at read 38. Diamonds are the ramp after droop 
subtraction (the saturated reads are not shown since they are not 
used in subsequent processing steps). The quadratic fit (used to de- 
rive the electronic nonlinearity correction) to these points is shown 
with the solid line. The final, linearized ramp is shown with the 
open squares. 



We also carried out detailed numerical modeling of the 
behavior of the Ge arrays. This modeling allowed ef- 
fects found in the laboratory testing to be investigated 
in more detail. For example, the modeling was able to 
show that the difference in hook behaviors between the 
70 an d 160 um arrays wa s due to their different illumina- 
tions l)Haegel et al.ll200ll) . The numerical modeling was 
also crucial to the understanding of the behavior of small 
signals on the detectors. For example, this modeling was 
able to determine that the stim flash latents (see ^5.2. 1[) 
were additive, not multiplicative. This understanding 
then guided the efforts to remove this signal. 

5. RAMPS TO SLOPES 

The first part of the processing fits the ramps to pro- 
duce slopes for each DCE. The processing for the Ge (70 
and 160 /im) and Si (24 /xrn) RAW mode data is similar, 
differing only in the instrumental signatures removed. 
First, reads which should be rejected from the linear fits 
are identified. Reads are rejected if they represent miss- 
ing data, autoreject reads, or saturated data. Second, 
the ramps are corrected for instrumental effects. These 
are dark current (Si only), rowdroop (Si only), droop 
(Si only), electronic nonlinearities, and stim flash latents 
(Ge only). Third, jumps in the ramps usually caused 
by cosmic rays are identified. In the process, reads that 
are abnormally noisy are identified as noise spikes. Fi- 
nally, all the continuous segments in each ramp are fit 
with lines and the resulting slopes averaged to produce 
the final slope for each pixel. An example of a 24 /im 
ramp is given in Fig. The 70 and 160 /xm ramps are 
very similar to the 24 /zm ramp, except they do not have 
droop. The graphical representation of the data process- 
ing shown in Fig. 01 gives the ordering of the reduction 
steps. The processing for the Si SUR mode data is nec- 
essarily different as the ramps are fit on-board and only 
the slope and first difference images are downlinked. The 
following subsections will describe the Si and Ge RAW 
mode processing followed by a description of the neces- 
sary differences for the Si SUR mode processing. 



5.1. Steps Common to Si & Ge RAW Modes 
5.1.1. Rejected Reads - Autoreject and Saturated Reads 

There are two reasons to automatically reject (autore- 
ject) reads; to avoid reset signatures and to not use the 
ramps beyond 2 MIPS seconds for stim flash DCEs. All 
MIPS arrays are reset at the beginning of a ramp, and 
this has been seen to leave a signature in the first few 
reads. In general, this reset signature only affects the 
first read. The first read is automatically rejected for all 
three arrays. This is even true for SUR data for which 
the line fit is done on-board Spitzer. The 70 and 160 /im 
arrays can be operated with a reset in the middle of the 
DCE to improve performance. When this mode is used, 
the reset signature has been seen to last for 4 reads and 
these 4 reads are then automatically rejected. In a stim 
flash DCE only the first 2 MIPS seconds of a ramp are 
valid. After 2 MIPS seconds, the stim is turned off and 
after 2.5 MIPS seconds, a reset is applied. 

Finally, all reads that are below or above the allowed 
limits for the MIPS analog-to-digital converters (ADC) 
(soft saturation) or saturating the 70 and 160 /im read- 
out circuits (hard saturation) are flagged as low or high 
saturation, respectively. 

5.1.2. Electronic Nonlinearity Correction 

All three MIPS arrays display nonlinearities that have 
been traced to the electronics. For the 24 array 
these nonlinearities are mainly due to a gradual debi- 
asing which occurs as charge accumulates in each pixel 
during an exposure. For the 70 and 160 fim arrays, the 
readout circuits have been constructed to keep the same 
bias voltage across the detectors even as charge accumu- 
lates. Nevertheless, electronic nonlinearities arise due to 
the simplified CTIA circuit. 

The behavior of the electronic nonlinearities was deter- 
mined from extensive ground-based testing on the flight 
arrays. For the 24 firm array, the functional form was 
characterized from RAW mode data ramps; a typical case 
is shown in Figure [S] The ramps for most of the pixels 
can be nearly perfectly described by quadratic polyno- 
mial fits; the linear component of the fit gives directly the 
linearized signal. For the 70 and 160 /im arrays, the elec- 
tronic nonlinearities have been shown generally to have 
a quadratic shape with significant deviations. Correc- 
tions were tabulated as a lookup table to allow for the 
semi-arbitrary forms. For the 24, 70, and 160 /im arrays, 
the maximum nonlinearity at full well (ADC saturation) 
ranges over the array from ~ 10 — 15%, ~ 1 — 2%, and 
~ 0.5 — 1%, respectively. 

5.1.3. Ramp Jumps - Cosmic Rays, Readout Jumps, 
and Noise Spikes 

The main reason discontinuities or jumps appear in 
MIPS ramps is cosmic rays. Cosmic rays strike the Ge 
detectors (70 and 160 fim arrays) at a rate of one per 
pixel per twelve seconds. The rate on the Si detector 
(24 /im array) is much lower, due to its smaller pixels. It 
is also possible to get a ramp jump due to an anomaly we 
have termed a readout jump. Ground-based testing has 
shown that the entire output of one of the 32 readouts 
(4x8 pixels) on the 70 /im array occasionally jumps 
up and then jumps back down by the same DN amount 
approximately 1 second later. 
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Jumps in the ramps are detected using a combination 
of two methods. First, (n — 1) 2-point differences are 
constructed from the n reads and outliers are flagged as 
potential ramp jumps using an iterative sigma clipping 
algorithm. These potential jumps are tested to see if they 
are noise spikes or actual ramp jumps by fitting lines to 
the segments on either side of the potential jump. If the 
two fitted lines imply a jump that is smaller than the 
expected noise, then the jump is actually a noise spike, 
not a cosmic ray or readout jump. 

Secon d, a more sensitive tes t for ramp jumps is per- 
formed ijHesselroth et all 1200(1) . This method works by 
assuming each read in a ramp segment has a ramp jump 
after it and fitting lines to the resulting two subsegments 
on either side. The most significant ramp jump in the 
segment implied from the two line fits is then tested to 
see if it is larger than the noise. If so, then this read is 
labeled as a ramp jump. The process can be repeated on 
the subsequent ramp segments until no more jumps are 
found or a preset number of iterations have been per- 
formed. As this second method is more sensitive than 
the first, but significantly more computationally inten- 
sive, we combine the two methods to achieve the best 
sensitivity to ramp jumps with the least computation 
time. 

We explored the signatures of cosmic rays in ramps us- 
ing several hours' worth of 70 and 160 /im array data that 
were subject to constant illumination. We then extracted 
those ramps where we detected ramp jumps (assumed to 
be due to energetic particle impacts) and assessed the ef- 
fects on the ramp after the impact. On the 70 /jm array 
we find two main effects: a steepening of the ramp that 
lasts for a few reads and a persistent responsivity increase 
of ~ 1% after a big hit (see Fig. EJ. This is consistent 
with the slow responsivity increase observed during the 
radiation run. These results dictate our strategy for deal- 
ing with cosmic ray hits on this array: several reads after 
a hit should be rejected from slope fitting to ensure that 
the fast transient does not bias the slope measurement, 
while the small responsivity increase after large hits will 
be tracked by the stim flash measurements. 

The 160 /im array response to cosmic rays is somewhat 
different. We detected no fast transient within the ramp, 
but the slope of the ramp after a hit was often different 
from the slope before the hit. This slope change typi- 
cally did not persist into the next DCE, after a reset had 
occurred, as shown in Figure El Thus, we were unable 
to detect a persistent responsivity increase due to par- 
ticle impacts, in contrast to the accelerator data f H3.3() . 
Given that we are unable to predict how the slope will 
change after a particle hit and that the slope returns to 
its previous value after the next reset (usually the next 
DCE), the conservative strategy for dealing with particle 
impacts on this array is to simply ignore all data between 
a particle hit and the next reset. 

5.1.4. Line Fitting 

Slopes are determined for each ramp by fitting lines to 
all the good segments in a ramp. The slope for a ramp is 
then the weighted average of the slopes of the ramp seg- 
ments. The weight of each segment is determined from 
the uncertainty in the segment slope as discussed in the 
next paragraph. Each good segment of a ramp is iden- 
tified as containing only good reads and not containing 



any ramp jumps. Lines are fit to these segments with 
the standard linear regression algorithm. 

Calculating the uncertainties on the fitted slope and 
zero point is not as straightforward. The uncertainties 
on each read have both a correlated and random compo- 
nent. The correlated component is due to photon noise, 
as the reads are a running sum of the total number of 
photons detected. The random component is the read 
noise. We have derived equations for the linear fit uncer- 
ta inties for the c orrelated component following the work 
of lSparka (1998). The details of this derivation are given 
in appendix ^ The slope and zero point uncertainties 
are calculated for the correlated and random read uncer- 
tainties separately and then combined in quadrature to 
get the final uncertainties. 

5.2. Steps for Ge Raw Mode Only 
5.2.1. Stim Flash Latent Correction 

The calibration of the 70 and 160 /im arrays is di- 
rectly tied to the stim flashes measured approximately 
every two minutes. The brightness of these stim flashes 
is set as high as possible to ensure the best calibration 
(cf. i !3.2f> . These stim flashes produce a memory effect, 
called a stim flash latent, that is persistent for a brief 
time. Intensive measurements of stim flash latents have 
been performed at the University of Arizona on the 70 
and 160 /im characterization arrays. We determined the 
time constants, amplitudes, variations with the back- 
ground, and repeatability of the stim flash latents as well 
as the accuracy of the correction and the effects on the 
calibration of sources observed during the latent. 

To characterize the decay behavior of the latents, we 
fit an exponential law to the time signal of each array 
pixel. Each cycle is divided by the stim amplitude value, 
to have dimensionless data (fraction signal/stim). The 
function F used to fit the latent is a double exponential: 

F{t) = b + ai e~ t/Tl - a 2 e~ t/T2 (2) 

where t is the time after the stim is turned off, b is the 
background level, a\ and a 2 give the component ampli- 
tudes, and t\ and t 2 give the time constants. 

At 70 /im, only a single exponential is needed (thus 
a 2 = 0). The amplitude ai is always less than 3% of the 
stim amplitude, and in most of the cases below 0.5%. 
The time constant t\ ranges from 5s to 20s. As a func- 
tion of increasing background, a\ increases and t\ de- 
creases. The latents are repeatable to 15% or better. An 
example of the stim latent of one pixel on the 70 /im 
characterization array is given in Fig. [7] 

At 160 /im, the latency effect is more pronounced than 
at 70 /im. The amplitude a\ is less than 5% of the stim 
amplitude. The time constant t\ ranges from 5 to 20s. 
The amplitude a 2 is less than 3%. The time constant 
t 2 equals 20s at high background, and is negligible at 
low background. The amplitude a\ and time constant 
n are almost insensitive to the background. The latents 
are repeatable to 20% or better. Fig. also gives an 
example of the stim latent of one pixel on the 160 /im 
characterization array. 

In general, the stim flash latents are negligible ~30 sec- 
onds after the stim is turned off. In the first 30 seconds, 
the calibration of a point source might be overestimated 
by 1% at 70 /im and 12% at 160 /im if no correction is ap- 
plied. To correct for the stim latent contribution to the 




Fig. 6. — Each panel shows the ratios of the slopes of ramp segments within a DCE in which a particle hit occurred (designated la before 
the hit and lb after the hit) against the ratios of the slopes of ramps in the DCEs before (designated DCE 0) and after (designated DCE 
2) the DCE with the hit. The arrays were subject to a constant illumination level and so the points should be clustered around 1 if there 
was no effect due to particle hits. The 70 fim array is shown in panel (a) while the 160 [ira array is shown in panel (b). 
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Fig. 7. — Examples of stim flash latents in the 70 fim (left) and 160 fim (right) characterization arrays. These plots represent time series 
of a pixel output averaged over 50 cycles. Each 10s DCE has been subdivided in 5 sub-DCEs of 2s. The diamonds give the pixel signal 
during the stim DCE and the crosses the pixel signal during science DCEs. The solid line gives the fit. The 70 /im fit parameters for pixel 
(2,22) are: stim flash signal = 18004 DN/s, b = 372 DN/s (2.1% of stim flash), oi = 256 DN/s (1.4 % of stim flash), and n = 14s. The 
160 fim fit parameters for pixel (9,0) are stim flash signal = 26985 DN/s, b = 193 DN/s (0.7% of stim flash), ai = 1091 DN/s (4.0 % of 
stim flash), and t\ = 4.3s. 



pixel signal, we apply a time-dependent correction at the 
ramp level. We subtract the latent contribution, which 
is obtained by integrating Eq. O On pre-flight data, the 
amplitude of the latents after correction is reduced by a 
factor of ~2 at 70 /xm and ~4 at 160 ^m. 

5.3. Steps for Si RAW Mode Only 

5.3.1. Rowdroop Subtraction 

The rowdroop effect manifests itself as an additive con- 
stant to each individual pixel and is proportional to the 
sum of the number of counts measured by all pixels on 
its row, where a row is in the cross-readout direction. 
This effect is not completely understood, and is simi- 
lar to (but separate from) the droop phenomenon (see 
H5.3.2|) . The additive signal imparted to each pixel on a 
row is constant and exhibits no gradient or dependence 



with pixel position, thus, it is not related to a charge- 
bleed, or "Mux-bleed" effect. The rowdroop contributes 
a small amount to the flux of an individual pixel, and 
will only significantly affect pixels on rows with high- 
intensity sources. An example of rowdroop from ground- 
based testing is shown in Fig. |S] 

Using images of pinhole sources obtained in ground 
testing, we have computed the row droop constant of 
proportionality, A" r( j. This is the factor that gives the 
fraction of the total counts in a row which is the result 
of row droop and should be subtracted from each pixel 
in that row. We find that the constant of proportionality 
for the MIPS 24 array is K rd = 7.6±2.5x 10" 5 . Thus, 
the rowdroop contributes « 1% of the total number of 
counts on a row. The rowdroop is corrected for on a 
read-by-read basis. 
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Fig. 8. — Example of the rowdroop effect on the 24 fim array. Two pinhole sources were illuminated by an external stimulator source, 
producing point sources on the array. On the left panel, the rowdroop signal produced by the pinhole sources can be seen as a horizontal 
stripe of artificial signal (other structure around the point sources, such as the arc around the right-hand source, is due to reflected light 
from the pinhole apparatus, while the other apparent point sources on the upper half of the array are latent signals produced by earlier 
exposures of the pinhole sources at different positions). The right panel shows the same image, after the rowdroop correction has been 
applied. 



5.3.2. Droop Subtraction 

Droop is a constant signal added to each pixel by the 
readouts. The exact cause of droop is unclear. This ex- 
traneous signal, akin to a DC offset, is directly propor- 
tional to the total number of counts on the entire array at 
any given time. We have measured the constant of pro- 
portionality from ground test data. The droop coupling 
constant was measured to be 0.33 ± 0.01, which agrees 
well with the 0.32 determined by BNA. 

The droop correction algorithm first computes the 
mean signal on the array, which is then multiplied by 
the droop coupling constant to derive the droop signal, 
as given by 

12 ij Cd 



F,= 



i + cv 



(3) 



where Fd is the droop signal, is the signal on each 
pixel (comprising both the actual incident flux and the 
droop), N pix is the number of pixels, and Cd is the droop 
coupling constant. The resultant droop signal is then 
subtracted from the original signal on each pixel. 

Under normal circumstances, the uncertainty associ- 
ated with this process is at the ~ 1% level, limited mainly 
by the uncertainty on the coupling constant. However, 
greater uncertainties arise when pixels are saturated; 
since ADC saturation occurs well before hard detector 
saturation, droop signal will still accumulate for an inci- 
dent flux above the ADC saturation level. In this case, 
the actual signal ramp must be extrapolated beyond the 
saturation point. The droop signal is determined by ex- 
trapolating a fit to the unsaturated portion of the ramp. 
As with the rowdroop correction, the droop correction is 
done on a read-by-read basis for RAW mode 24 /j,m data. 

5.3.3. Dark Subtraction 



Dark subtraction is done at each read using a dark 
calibration image containing the full dark ramp for each 
pixel. This step serves both to remove the (small) dark 
current contribution and the offset ramp starting points, 
so that each ramp starts near zero. 

5.4. Steps for Si SUR Mode Only 

The majority of the 24 /im data are taken in the SUR 
mode instead of the RAW mode. In the SUR mode, a 
line is fit to the data ramp on-board the spacecraft. The 
resulting slope and first difference (difference between the 
first two reads of the data ramp) images are downlinked 
instead of the full ramp. The first difference frame effec- 
tively increases the dynamic range of the SUR mode as 
signals that saturate somewhere in the ramp, but after 
the second read, will have a valid measurement in the 
first difference frame. To reduce the data downlinked, 
any first difference value that is from a ramp which does 
not saturate is set to zero. This increases the compress- 
ibility of the first difference frame. 

5.4.1. SUR Saturation Detection 

There can be degeneracy of SUR slope values due to 
the possibility of saturation. The possible slope value 
for a given pixel reaches a maximum at full well, the 
point of ADC saturation. After that point, as the data 
ramp reaches saturation at the last few reads, the slope 
value will begin to decrease because the on-board SUR 
algorithm does not reject saturated reads. In cases of 
extreme saturation, the slope becomes quite small, and 
can eventually become zero if saturation occurs within 
the first few reads. The first difference value is provided 
to break this degeneracy. We have employed a conserva- 
tive threshold value for the first difference, above which 
a pixel is flagged as being likely saturated. ADC satura- 
tion occurs at +32768 DN (see Figure [5] for an example 
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of a saturated RAW ramp) . Assuming a linear ramp, the 
first difference for a ramp that just saturates on the final 
read would be 65536/n rea£ 2, where n rea d is the total num- 
ber of reads in the data ramp. For example, there are 60 
reads in a 30 second DCE, yielding an ideal saturation 
threshold of ~ 1100 DN/read. To be more conservative, 
we actually employ a threshold value of 1000 DN/read 
for a 30 second exposure time and scale this for other 
exposure times. Since the data ramps arc not linear, the 
actual first difference threshold is larger than our chosen 
default value, so most cases of saturation will be flagged. 
The only exception being saturation at the first read, in 
which case both the slope and the first difference would 
be zero. For all pixels that have been flagged for satura- 
tion, the first difference value should be used in place of 
the slope. 

5.4.2. Rowdroop, Droop, and Dark Subtraction 

The rowdroop and droop subtraction is done in the 
same way for SUR mode as for RAW mode, except that 
the corrections are performed on the slope and first dif- 
ference images. 

5.4.3. Electronic Nonlinearity Correction 

Because the SUR data do not preserve the actual data 
ramps, the linearity correction made somewhat compli- 
cated. Nevertheless, the quadratic behavior of the ramps 
can be used to analytically determine the linearization 
of the SUR slope values. This correction depends on the 
observed SUR slope value, exposure time, and known 
quadratic nonlinearity. Note that saturation invalidates 
this method, as the SUR slope-fitting algorithm does not 
reject saturated reads. In this case, no linearity correc- 
tion is applied. 

6. SLOPE IMAGE CALIBRATION 

The next step in the MIPS data reduction is to cali- 
brate the slope images while removing instrumental ef- 
fects with time constants longer than the DCE expo- 
sure times. The instrumental effects corrected at this 
stage include latents (24 /xm), responsivity drift (70 and 
160 /im), pixel-to-pixel responsivity variations, the tele- 
scope illumination pattern, and flux nonlinearities (70 
and 160 /im). The graphical representation of the data 
processing shown in Fig. 0] gives the ordering of the re- 
duction steps. Since the time dependent responsivity of 
the Ge arrays requires additional calibration steps than 
is usual for more common array detectors, we give the 
mathematical basis of our Ge slope calibration in iJ6.ll 

6.1. Principles 

Ignoring the 70 and 160 /im flux nonlinearities, an un- 
calibrated slope image can be represented by 

U(i,j,t n ) = [I(i,j)0(i,j)+D(i,j)]R(i,j,t n ) (4) 

where I(i,j) is the science image of interest, 0(i,j) rep- 
resents the telescope and instrument optics (the mean of 
0{i,j) is one), D(i,j) is the dark current, and R(i,j,t) 
is the instantaneous responsivity of the array; i,j rep- 
resent the pixel coordinates and t n the time of the 
nth DCE. Calibration involves isolating the flux 

from the sky+object in the above equation. The term 
0(i, j)R(i, j,t n ) is the equivalent of a traditional flat- 
field term. As R(i,j,t n ) is a rather sensitive function of 



time for the 70 and 160 /im detectors, however, a global 
"flat-field" cannot be determined, but must be derived 
for each DCE separately. The stimulators provide the 
means to monitor R(i,j, t n ) and all science observations 
will be bracketed by stim flashes. Stim flash images will 
be equivalent to science frames with the addition of a 
stimulator illumination pattern: 

U s tim,N — [S(i, j)0(i,j)+D(i,j)]R(i,j, t N ) (5) 

where S(i,j) is the illumination pattern introduced on 
the array by the stim flash with the mean of S(i,j) equal 
to one. MIPS observations include the requirement that 
each stimulator flash will be preceded by a background 
exposure with the identical telescope pointing; thus for 
the N th stimulator DCE there exists a background DCE 
taken at time tjy — e, 

U bkgd , N = [I(i,j)0(i,j)+D(iJ)]R(iJ,t N -e). (6) 

If we assume that the responsivity of the array R(i,j, t) 
doesn't change dramatically between times t^ and tj^ — e, 
i.e. R(i,j,ti\r) ~ R(i, J, ijv — e), we can construct for each 
stimulator flash a background subtracted stim flash: 

Ustim.N - Ubkgd,N-e ~ S(i, j)R(i, j,t N ). (7) 

With background subtracted stim flashes determined 
from Eq. for all stim flashes in the data set, an in- 
stantaneous stim can be determined for any time, t n , by 
interpolation from bracketing stim flashes: 

S(i,j)R(i,j,t n ) = F[S(i,j)R(i,j,t N )] (8) 

where F[] is some interpolating function on background 
subtracted stims for times tN bracketing t n . Analysis of 
Ge characterization array data indicates that a weighted 
linear fit (weighted by the uncertainty in the stim flash 
frames) to two stim flashes on either side of the data 
frame (a total of four stim flashes) provides the optimal 
strategy for determining the instantaneous stim ampli- 
tude (repeatability to ~1% on most backgrounds). Di- 
viding science frames, Eq. by the interpolated instan- 
taneous stim, Eq. [SJ produces 

U data (i,j) = [I d ata(i,j)0(ij) + D(i,j)]/S(i,j). (9) 

While we have removed the time dependent responsivity 
variation, the data of interest, Idata{i,j), are still mod- 
ified by the optical response and the dark current; in 
addition, we have introduced the stimulator illumination 
pattern into our data. Fortunately, since the time de- 
pendence has been removed, we can remove these other 
instrumental signatures through carefully accumulated 
calibration data. 

First, the dark correction, D(i,j), can be determined 
from a sequence of exposures as above, with the addi- 
tional constraint that the scan mirror be positioned such 
that no light from the "sky" falls on the detector. Thus 
the data and stim flashes in a dark current data sequence 
are represented by 

U(i,j,t n )=D(i,j)R(i,j,t n ) (10) 

and 

U stim = [S(i,j)+D(i,j)]R(i,j,t N ), (11) 

respectively. The dark data are corrected for responsiv- 
ity variations exactly as described above and the indi- 
vidual frames combined to produce an average dark cur- 
rent, D(i,j)/S(i,j). Subtracting this dark current from 
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science frames that have been corrected for responsivity 
variations, Eq. yields 



Udataihj) 



Idata{i,j)Q(i,j) 

S(i,j) 



(12) 



our responsivity, dark corrected science frame. What 
remains is to correct for the telescope optics, 0(i,j), and 
the stim illumination pattern, S(i,j). 

Correcting for the combined illumination pattern of the 
telescope and stim involves a standard series of MIPS ex- 
posures, i.e. data frames interspersed with stim flashes. 
As such, they may be represented by equations of the 
form Eq. where the I n (i,j) represent dithered images 
of "blank" sky fields. Calibrating the sequence by cor- 
recting for responsivity variations and dark current as 
above results in a series of images 



In(i,j)Q(i,j) 
S(i,j) ■ 



(13) 



Since by construction, the I n (i,j) are dithered images 
of "smooth" regions, if a large number of I n (i,j) are 
acquired, they may be median combined to remove point 
sources (and extended sources if dithered "sufficiently"), 
cosmic rays, etc. resulting in 

U iUum {i,j) = (I n (i, J)>SM = (14) 

where 0(i,j) and S(i,j) are constant regardless of tele- 
scope pointing. Hence the median only affects the chang- 
ing sky image as the telescope is dithered. The constant 
C in equation 1151 may be set to one resulting in the illu- 
mination correction frame 



Q(i,j) 
S(i,j) 



(15) 



The responsivity corrected, dark subtracted data (equa- 
tion I12f> are now divided by the illumination correction 
resulting in 



Udata(i,j) 



I(i,j)Q(i,j) , 0(1,3) _ r( . 

s(ij) 's(i,j) - 1{h3h 



(16) 



and we have recovered the quantity of interest, the as- 
tronomical sky I(i,j). Suitable observations of standard 
stars can then be used to convert instrumental counts to 
physical units (e.g. Janskies). 

6.2. Dark, Flat Field, and Illumination Correction 

The dark, flat field, and illumination correction cali- 
bration images described above will be obtained through- 
out the life of the mission. Example prefiight calibra- 
tion images are shown in Figs. 1911 II Simulations indi- 
cate that high S/N flat field and illumination correction 
images (~ 0.5% RMS) can be obtained with dithered ob- 
servations of "smooth" areas of the sky: ^60-100 DCEs 
at 24 //m and ~200 DCEs at 70 /itm are required. At 
160 /im, the situation is less ideal, with simulations indi- 
cating as many as 500 DCEs may be required to produce 
flats to better than 1% RMS. 

6.3. Si Latent Correction 

Si IBC arrays are known to have considerable latency, 
where the signal induced by bright illumination persists 



after the illumination has terminated. Ideally, if one 
knows the position of a source exposed on the array and 
the latency decay behavior, these artifacts can be sub- 
tracted from an image. We have characterized the latent 
behavior from ground test data. Several different con- 
ditions were explored, including varying brightnesses of 
the illuminating source, varying brightnesses of the back- 
ground, initial bias boosts, and changing the number of 
resets via different exposure times. A bias boost can 
flush out most of the trapped charge, but resets are not 
nearly as effective. Since bias boosts will only be done in 
the first DCE of each observation, we correct for latent 
residuals in the data processing. 

The latent decay curve can be described by single ex- 
ponential, given by 



m(t) 



m 



pe 



-t/r 



(17) 



where m is the slope in the absence of a latent, p is the 
initial value of the latent, and r is the latent time con- 
stant. Based on the limited ground data, the latent pa- 
rameters (p and t) appear to be functions of background 
levels, number of resets and possibly location on the ar- 
ray. In general, the latent contribution is about ~ 1% of 
the initial source brightness ~5 sec after that source has 
shut off. Higher background yield slightly higher values 
for p and lower values of r. The value of r is in the range 
of 12 ± 5 seconds. 

6.4. Ge Flux Nonlinearity Correction 

Both the 70 and 160 /im arrays exhibit nonlinearities 
that are dependent on the incident point source flux as 
well as the background. These are termed flux nonlin- 
earities and have been observed in data taken with the 
characterization array as well as the flight array. As is 
usual for the Ge arrays, each pixel shows flux nonlin- 
earities with a different dependence on source flux and 
background. Correcting for this effect can be broken into 
two pieces: (1) removing the pixel to pixel differences 
in the nonlinearity followed by (2) the application of a 
global nonlinearity correction as a function of the source 
brightness and background. 

The pixel to pixel variations in the flux non-linearity 
may be mapped by analyzing the ratio of two stim 
flashes, where one is the standard on orbit calibrating 
stim flash. Measured differences in the ratio from pixel 
to pixel can be used to correct each pixel to the same 
flux non linearity for the given background and source 
(second stim flash) amplitude. Repeating the measure- 
ment for a variety of second stim flash amplitudes (up to 
saturation for each pixel) and backgrounds will map out 
the correction. The second, global, stage of the correc- 
tion can be characterized by observations of calibration 
stars with a range of known brightness ratios on similar 
backgrounds. The combination of these two tasks out- 
lined above should provide a good measurement of the 
flux nonlinearity correction for a range of backgrounds. 
This correction will improve continuously during the mis- 
sion as the range of backgrounds and calibration stars 
expands. 

6.5. Flux Calibration 

The absolute calibration of MIPS will rely on a 
well determined anchor at 10.6 /im using the fun- 
damental calibrators a Boo, a Tau, and (3 Gem 
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Fig. 9. — Calibration Images for the 24 fim array. They are the dark (left) and fiat field (right). The displayed range for the dark image 
is to 4 DN/s. The displayed range for the flat field image is 0.9 to 1.1. The pixels displayed as solid black correspond to the handful of 
known bad pixels. 




Fig. 10. — Calibration Images for the 70 ^m array. They are the dark (left), the illumination correction (middle), and an example of a 
stim flash (right). The displayed range for the dark image is to 0.02. The displayed range for the illumination correction image is to 
2.0. The displayed range for the stim flash is to 40,000 DN. 




Fig. 11. — Calibration Images for the 160 (im array. They are the 
dark (top), the illumination correction (middle), and an example 
of a stim flash (bottom). The displayed range for the dark image 
is to 0.65. The displayed range for the illumination correction 
image is to 2.0. The displayed range for the stim flash is 4,000 
to 24,000 DN. 



URieke. Lebofskv. fc Low! Il985t ICohen et all llflfl"^ . 

Three independent methods will be used to extrapolate 
the calibration at 10.6/im to the MIPS bands: (1) 
Solar analogs, (2) A star atmospheric models, and (3) 
semi-empirical models of K giants. Grids of stars for 
each method have been observed from the ground and 



tied to the fundamental calibrators at 10.6 /im. For 
the solar analog stars, on orbit observations at 24, 70, 
and 160 ^m are being compared with extrapolations of 
empirical measurements of the sun extrapolated into the 
MIPS bands. A grid of A stars has been observed in all 
three bands on orbit and compared to extrapolations of 
A star atmosphere models to the MIPS bands. While 
the solar analog and A star calibrators will be observed 
in the MIPS 160 /im band, the K giant calibrators will be 
the only ones detectable at high signal-to-noise in that 
band. On orbit observations of the K giant calibrators 
are being compared to theoretical e x trapo l ations o f 
model atmospheres, eg. ICohen et alJ 1119951 H996albD 
extrapola ted to longer wavelengths using the Engelke 
function (Engelke 1992). Absolute flux calibrators will 
be observed throughout the lifetime of the mission. 

7. USING REDUNDANCY TO IMPROVE CALIBRATION 

The last step in the reduction of MIPS data is to use 
the redundancy inherent in the observing modes to im- 
prove the removal of instrumental signatures. This step 
is mainly for the 70 and 160 /im data due to the chal- 
lenging aspects of Ge detector calibration. We define the 
level of redundancy to be the number of different pixels 
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that measure the same point on the sky. Our approach 
will be to look for known instrumental signatures (as a 
function of time) in the difference between what a partic- 
ular pixel detects and what all the other pixels detected 
for the same sky locations. This is possible because the 
observing strategy has been designed so that each point 
on the sky will be observed multiple times by different 
pixels. 

Table ^ shows the minimum level of redundancy for 
each MIPS observing mode. Many MIPS observations 
are taken with multiple cycles resulting in significantly 
higher redundancies. It is recommended to have a mini- 
mum redundancy of four. 

7.1. Algorithm 

The basic algorithm for using redundancy to refine the 
instrumental signature removal is as follows. 

1. Create a mosaic of all the images in question. Dur- 
ing the mosaic creation, use a sigma rejection al- 
gorithm to remove data that are deviant from the 
majority of the observations. 

2. Use the mosaic as a "truth" image of what each 
image should have measured. 

3. For each pixel, subtract the actual from the "truth" 
measurements to create a measurement of the time 
history differences. 

4. Examine the difference time history for known in- 
strumental signatures. While many instrumental 
signatures could be present, we plan to concentrate 
on stim latent residuals and systematic differences 
between "extended" sources and point sources. Ac- 
tual on-orbit data will guide the details and number 
of instrumental signatures that are corrected using 
redundancy. 

5. Correct for all instrumental signatures that are 
found to be significant. 

6. Iterate steps 1-5 until no new significant instrumen- 
tal signatures are found. 

The input to this algorithm is calibrated slope images. 
The output product of this algorithm is enhanced images. 
A useful side product will be the mosaicked image of the 
object. 

7.2. Distortions of Arrays 

To use the redundancy to remove additional instru- 
mental signatures we must first coadd all related obser- 
vations into a single mosaic. Because the MIPS optical 
train is made up of purely off-axis reflective elements 
there exist scale changes and rotations across the re- 
imaged focal plane. To coadd images taken at different 
places on the array, it is crucial to correct the data for 
these distortions. 

We used the Code V optical models for Spitzer/MIPS 
to estimate the distortions present in the images from the 
three MIPS detectors. The results from Code V allow us 
to determine distortion polynomials which can then be 
used to correct for the distortions. We estimated the dis- 
tortions by setting up a grid of equally spaced points in 



i 



1 1 1 

, M .1 ! 

MtttUl ill 
\\\\ \ ' 

\\\\\"" 

\\ \ \ yv-i 

x \ v \ -i ^ 
w v v \ v ^ 
x\ % v 

x> V v \ 



1,0 pi*el 



_1_ 



* * v \ 
v v ^ \ V 
^ \ \ v v 

■ r l i \ \ r \ \ v 

f 1 -T f \ t \.<\ \ \ 

r 1 1 I I 1 \ \ \ \ 

r i i i l l l \ \ \ 

I''' 



20 

Cross Scan 



Fig. 12. — 70 fim narrow field mode residuals between the dis- 
torted and undistorted points after the object angles were con- 
verted to pixels. Note the 1 pixel scale in the lower left-hand cor- 



the field of view at a specific scan mirror angle. The chief 
ray from each object point was traced through the sys- 
tem to where it was imaged on the focal plane. In a per- 
fect optical system the image points would map perfectly 
from the object with a possible change in magnification. 
The difference between the ideal location and the actual 
location is the distortion. For example, Figure El is a 
vector plot of the distortions present in the 70 /im nar- 
row field array. The equally spaced grid of points present 
the focal plane points and the ends of the vectors corre- 
spond to the object points, after a plate scale factor was 
applied. The difference in the points (the length of the 
vector) is caused by the distortions. 

Table lists the scale change of the field of view of 
the different MIPS arrays. The scale change is defined 
as (maximum length of distorted field - minimum length 
of the distorted field) /(minimum length of the distorted 
field). From a distortion standpoint, it is useful to look 
closely at individual pixels to see how distortion changes 
the area imaged on the pixel. Figure ED is a plot of a 
distorted pixel in the 70 //m narrow field array. One 
can see that the distorted pixel changes shape from a 
square to a somewhat trapezoidal shape. The ratio of 
the distorted to undistorted pixel area is 1.19. Table 01 
lists information on how distortion affects the area im- 
aged on individual pixels on the different arrays. The 
distorted pixel area ratio is defined as (distorted pixel 
area) /(undistorted pixel area). 

Following the procedure of converting the pixel 

coordinates to wo rld coordinates outlined in 

iCreisen fc Calabrettal (|2002t) . the distortion correc- 
tion is applied to the pixel coordinates before any 
other transformations. The distortion correction is 
accounted for by distortion polynomials. The distortion 
polynomials give the additive correction to map the 
distorted pixel coordinates, u,v to the distortion cor- 
rected pixel coordinates p, q. Thus, p = u + F(u, v) and 
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TABLE 1 

Redundancy in MIPS Observing Modes 



Type of AOR 


24 fim 


70 /im 


160 fj,m 


Photometry, compact 


14 


10 


2 


Photometry, large 


10 


6 


1 


Photometry, compact super resolution 


11 


8 


6 


Photometry, large super resolution 


10 


8 




Scan Map, slow 


10 


10 


1 


Scan Map, medium 


10 


10 


1 


Scan Map, fast 


5 


5 


0.5 


SED 




2 




Total Power 


1 


1 


1 



TABLE 2 
Distortions in MIPS 
Detectors 



Detector 


% FOV Scale 


24 fim 


2.84 


70 fim Wide 


0.2 


70 fxm Narrow 


7.70 


160 /im 


7.78 



32.75 



32.00 - 



„ 31.25 - 



30.50 - 



29.75 I i i I i i I i i I i i 

-0.50 0.25 1.00 1.75 2.50 

Cross Scon 

Fig. 13. — 70 /im NF mode distorted pixel. This pixel is located 
in the right hand corner of the array with (area distorted/area 
undistorted) = 1.1929. In this plot the distorted pixel is plotted 
with a solid line and the undistorted pixel is plotted with a dot- 
dash line. 



q = v + G(u, v), where 

F(u,v) = A2ou 2 +Ao2V 2 +Aiiuv+A 3 ou 3 +A2iu 2 v+A 1 2vu 2 + 

(18) 

and 

G(u,v) = B 20 u 2 +B () 2V 2 +B 11 uv+B 3 ( ) u 3 +B2 lU 2 V+Bi2VU 2 

(19) 

7.3. Mosaicking Details 

The ability to remove additional instrumental signa- 
tures is dependent on creating a high resolution mo- 
saicked image. For the mosaicked image to be of sufficient 



resolution, the mosaicked pixel sizes must be smaller than 
the original input pixels. While many mosaicking pro- 
grams compensate for undersampling by making use of 
dithered data, the MIPS data are well sampled and do 
not require this compensation. Instead, our focus is co- 
adding the related calibrated images into a single image 
without interpolating between pixels. Therefore, we al- 
ways work on the coordinates of the corners of a pixel, 
transforming them from the input image coordinate sys- 
tem to the output mosaicked coordinate system. The 
output mosaic image is on a single tangent plane. In the 
transformation of the pixel corners in the input image 
pixels to their location on the output mosaic image the 
corners are corrected for distortion, converted to right 
ascension and declination, and then projected onto the 
tangent plane defined by the right ascension and decli- 
nation of the mosaic center. Figure ITT1 is an example of 
three images which overlap each other on the mosaicked 
plane. In the process of establishing the location of the 
image pixel corners on the mosaicked plane, the link and 
overlap coverage between the input pixel and the out- 
put pixels it falls on is determined. A critical step in 
removing residual instrumental signatures based on the 
co- added mosaic image depends on correctly linking each 
mosaic pixel with each image pixel that overlaps it (and 
vice versa) and accurately determining the degree of over- 
lap. Essentially each output sub-sampled mosaic pixel 
becomes a cube of data, with each plane in this cube 
representing the information in each overlapping image 
pixel. The surface brightness and uncertainty associated 
with each mosaic pixel is found by weighted averaging 
the overlapping planes of data. In the surface brightness 
case, the weighting is based on the overlap coverage and 
uncertainty associated with the input image pixel. 

The information in each mosaic pixel is based on mul- 
tiple observations of a single area on the sky. This re- 
dundancy of data can be used to identify cosmic rays or 
any single image pixel measurement that deviates from 
thpexpected mean of multiple observations and expected 
r?oise. For example, for a 70 /im photometry observing 
mode cycle, if one of the pixels suffers from a much larger 
stim latent than the other observations, it will stand out 
B@$yi.be identified as an outlier. As an outlier, it will not 
be used in creating the mosaicked image. After all the 
outliers have been determined, then the links between the 
output mosaic pixel and the input image pixels are used 
to tally the number of times an image pixel was flagged 
as an outlier. If the majority of the time an image pixel 
was flagged as deviant, then this pixel is flagged in the 
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TABLE 3 
Distorted Pixel Area Ratio 







Area Rat 


io 




Detector 


Mean 


Std. Dev. 


Min. 


Max. 


24 /an 


0.9998 


0.0282 


0.9406 


1.0613 


70 /im Wide 


1.0027 


0.0042 


0.9973 


1.0148 


70 (im Narrow 


1.0129 


0.0664 


0.8913 


1.1929 


160 fim 


0.9781 


0.0361 


0.9007 


1.0137 



Image. 1 



Image. 2 



Image 3 







Mosaicked Grid 

Fig. 14. — An example of how a single pixel from three different 
images (only 4x4 pixels shown) are overlapped on the mosaic 
image is shown. The solid point represents the same location on 
the sky as it would be observed in each image. 



original data as an outlier. If a sufficiently large number 
(about 1%) of the input image pixels are flagged as out- 
liers then the mosaic step is repeated. The final output 
is a mosaic image that can then be used as the "truth" 
image of what each image should have measured. Fol- 
lowing the steps outlined in section 17. II this truth image 
is used to remove residual instrumental signatures. 

8. INITIAL TESTING WITH FLIGHT DATA 

With the successful launch of the Spitzer Space Tele- 
scope in August 2003, these reduction algorithms were 
tested against MIPS flight data of astronomical sources. 
This testing has validated the algorithms described in 
this paper, but has also shown that a number of mod- 
ifications will be needed to handle the realities of flight 
data. The i nitial results o f this testing are summarized 
here and in iGordon et alJ ( |2004f) , but a full accounting 
will be a subject of a future paper when the final MIPS 
reduction algorithms are known. 

There were two significant changes in the instrument 
operations which are not easily correctable by reduction 
algorithms. The 70 /im array was found to suffer from 
a cable short induced sometime between ground testing 
and flight. This short injects a large amount of noise into 
one half of the 70 /im array resulting in a useful array 
of only 16 x 32 pixels. The 160 /im array was found to 
suffer from a "blue-leak" caused by an unintended reflec- 



tion from the blocking filter which passes through the 
bandpass filter. This "blue-leak" results in an approxi- 
mately factor of 15 image leak for stellar sources. This 
leak means that asteroids are now the primary calibra- 
tors for 160 /im. Other than bright stars, the leak signal 
is below the confusion limit for most science targets as 
they have much smaller blue/160 /im ratios. 

At 24 /im, the pre-flight reduction algorithms were 
found to work well with only three changes needed. First, 
the Row Droop Correction does not seem necessary, but 
extensive testing has yet to be completed. Second, an 
additive offset in the second read of every ramp was 
found which produced a low level (1-2%) gradient in final 
mosaics. A straightforward correction for this has been 
implemented using RAW and SUR data for calibration. 
Third, scan mirror angle dependent flat fields are needed 
due to contamination of the scan mirror by small parti- 
cles. This contamination is seen as dark spots in indi- 
vidual images which move with scan mirror angle, but 
not spacecraft offsets. With these three modifications to 
the preflight algorithms, MIPS is producing high quality 
24 /im images which are well calibrated. 

At 70 /im, flight data has validated the basic structure 
of the preflight reduction algorithms but significant mod- 
ification is required to account for time dependent behav- 
iors. The stim flash latents were found to grow in am- 
plitude quickly after anneals. With a similar timescale, 
the residual background time dependence (after correc- 
tion using the stim flash amplitudes) was seen to grow. 
These two facts required hand reductions to remove the 
stim flash latents and background variations to produce 
good quality mosaics at 70 /im. These two effects are 
prime candidates for removal using redundancy, but the 
effectiveness of automatic removal has not been demon- 
strated yet. 

At 160 /im, the basic preflight algorithms have been 
validated from comparison with flight data. Some dif- 
ferences in detector behavior were seen in flight data. 
For example, the stim flash latents have a faster time 
constant than in preflight data. At this time, the nonlin- 
earities in the 70 and 160 /im arrays have not been well 
enough characterized with flight data to validate this sec- 
tion of the preflight algorithms. Finally, the cosmic ray 
rate seen in the Ge arrays has been seen to be about a 
factor of two over preflight predictions, 1 cosmic ray ev- 
ery 12 or so seconds. The ramp jump detection has been 
seen to work well and line segment fitting removes the 
majority of the effects of these cosmic rays. Some resid- 
ual effects remain and additional characterization may 
lead to algorithms to remove these additional effects. 

The effectiveness of the algorithms described in this pa- 
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Fig. 15.— Images of the 24, 70, and 160 (im observed PSFs are shown with dimensions of 60" X 60", 100" X 100", 150" X 150", 
respectively, for The 24 /an PSF was created from observations of HD 159330 a K2III star with a predicted flux of 0.570 Jy. The 70 fim 
PSF was created from fine-scale observations of HD 131873 a K4III with a predicted flux of 3.14 Jy. The 160 fim PSF was created from 
observations of Harmonia an asteroid with a predicted flux of ~1.5 Jy. 



per as well as the design of MIPS is attested by the point- 
spread-functions (PSFs) constructed from flight data at 
24, 70, and 160 fim shown in Fig. EJ These PSFs all 
clearly have a well-defined first Airy ring with the 24 fjxa 
PSF also exhibiting a well-defined second Airy ring. All 
three PSFs are well represen ted by the predictions of 
TinyTim models (|Krist»1993fl adapted to MIPS. In ad- 
dition, there are many papers written using MIPS flight 
data for the Spitzer Special Astrophysical Journal Sup- 
plement Issue (2004, ApJS, 154). 

9. SUMMARY 

This paper has described the preflight data reduction 
algorithms for all three arrays for the MIPS instrument 
on Spitzer. These algorithms have been guided by ex- 
tensive laboratory testing of the Si (24 /im) and Ge (70 
and 160 /im) arrays. In addition, numerical modeling of 
the Ge arrays has provided important insights into their 
behavior. 

The design and operation of the MIPS instrument has 
been summarized to give sufficient background for under- 
standing the data reduction algorithms. The design and 
operation of the MIPS instrument is mainly driven by 
the needs of the Ge arrays. As Ge detectors display sig- 
nificant responsivity drift over time due mainly to cosmic 
ray damage, the MIPS observing modes include frequent 
observations of an internal illumination source. In addi- 
tion, most MIPS operating modes have been designed to 
provide significant redundancy to increase the robustness 
of the MIPS observations against detector effects. 

The data reduction for the MIPS arrays is divided 
into three parts. The first part converts the data ramps 
into slope measurements and removes detector signatures 
with time constants less than approximately 10 seconds. 
These detector signatures at 24 /im include saturation, 
dark current, rowdroop, droop, electronic nonlinearities, 
and cosmic rays. At 70 and 160 /im, the detector signa- 
tures removed include saturation, electronic nonlineari- 
ties, stim flash latents, and cosmic rays. The resulting 
slopes are determined from linear fits and their uncer- 



tainties are computed accounting for both the random 
and correlated nature of the data ramp uncertainties. 

The second part of the MIPS data reduction converts 
the slopes to calibrated slopes and removes detector sig- 
natures with time constants larger than approximately 
10 seconds. At 24 /im, this translates to applying a flat 
field, correcting for object latents, and applying the flux 
calibration. At 70 and 160 /im, this step includes sub- 
tracting the dark, flat fielding using an instantaneous 
flat field, correcting for the flux nonlinearities, and ap- 
plying the flux calibration. A flat field specific to each 
70 and 160 /im image is required to correct for the 
time-dependent responsivity of the Ge arrays. It is con- 
structed from the frequent stim flashes and a previously 
determined illumination correction. 

The third data reduction step is to use the spatial re- 
dundancy inherit in the MIPS observing modes to im- 
prove the removal of instrumental signatures. This step 
is only applied to the Ge data. Known instrumental sig- 
natures are searched for in the difference between what 
a specific pixel and what all other pixels from the same 
sky locations detected. If instrument signatures are de- 
tected, they are removed and the process is repeated. 
This method is iterative in nature and will require care 
to avoid introducing spurious signals into the data. The 
design of this portion of the data reduction algorithms is 
necessarily the least developed because only after Spitzer 
launches will it be known which instrumental signatures 
are important to correct with this method. 

Finally, initial testing using flight data from MIPS 
has validated these data reduction algorithms, but some 
modification is necessary to account for the realities of 
flight. A future paper will describe these modifications 
in detail once they have been devised and tested. 
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APPENDIX 

LINEAR FIT TO DATA WITH CORRELATED AND RANDOM UNCERTAINTIES 



When a detector is non-destructively read out multiple times before resetting the resulting data ramps represent 
correlated measurements. This is because measurement yi+i is equal to yi +pi where pi is number of photons detected 
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in the time between yi and yi+i- This statement ignores the effects of read noise, which produces uncorrelated 
uncertainties on the yt measurements. While fitting lines to data with correlations is a complex subject, the form 
of the correlations in the case of non-destructive readouts allows analytic equations to be derived for the linear fit 
parameters and uncertainties. We present a derivation of equations for linear fit parameters and their uncertainties 
for the case of a data ramp with correlated reads and no read noise. This derivation is based on a similar derivation 
hv lSnaxksl fl99l for NICMOS data ramps but is slightly more general. As part of this derivation, it can be seen that 
the linear fit parameters derived assuming either random or correlated uncertainties are equivalent. This is not the 
case for the uncertainties on the fit parameters, which is the main motivati on for this derivation. 

The basics of fitting a line to data with random uncertainties are given in lBevington &: Robinson! l|1992f) . We repeat 
their results here, as the derivation for correlated uncertainties draws directly from this work. In fitting data to a line 
of the form 

y. k = a + bxi (Al) 

the fit parameters and their uncertainties are 



S XX Sy S x S x y 

A ' 

SSxy SxSy 



(A2) 



b= , (A3) 

a a {ran) 2 = ~^-, and (A4) 
a b {ran) 2 = ^ (A5) 



where ./V is the number of (xi,yi) measurements, cr(yi) is the uncertainty on each measurement of yi, 

N 1 

Y — 



N i 

S = Y,-7-T^ (A6) 



N 

5 ^E^i> ( A7 ) 

N 2 

S** = Y,-t%*> (A8) 



N 
N 

S - = E^ - d ( Al °) 

A = SS XX -(S X ) 2 . (All) 

These equations assume that the measurements of yi are independent. 

To determine the linear fit terms for a line fit to correlated data ramps, the assumption that the yi measurements 
are independent is not correct. The standard formulae need to be modified to sum over terms that are independent. 
The modifications start with realizing that 

Pi = Vi ~ Vi-i (A12) 
is the independent quantity in the absence of read noise. Any equation in the standard derivation that relies on the 
independence of yi needs to be modified to only depend on pi. Thus, 

JV 

s y=J2n% (Ai3) 



2/1 , 2/1 + P2 , 2/1 + P2 + P3 



(A14) 



ct(2/i) 2 ct(2/2) 2 cr(y 3 ) 2 

N 1 N 1 N 1 




(A16) 
(A17) 
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and using a similar derivation, 



S X y — 2_j 



Xi]Ji 



N 



N 



JV 



-111 



N 



i=2 
N 



k—i 



-yiS x + ^2 



i=1 



^{Vk? 



(A18) 
(A19) 
(A20) 



The standard equations IJA2I & IA3|) can then used to determine the best fit values of a and b for the case of correlated 
uncertainties. In fact, the values of a and b derived assuming correlated or uncorrelated uncertainties are exactly the 
same. The differences between the two types of uncertainties arises in determining a a and a b . 

To derive a n and ah for a data ramp with correlated measurements we start with equations 6.19 and 6.20 of 
iBevineton fc Robinsonl 1^)92). Converting from yi to pi as the independent variable gives 



N 



i=l 
N 

=£ 



dz 
dyi 



where z is either a or b. The partial derivatives needed are then 



da 
dpi 



1 



A 



dS. 



dpi 



dpi 



I ( N I N X 



and 



db_ 

dpi 



dS. 



as,. 



dpi x dpi 



Thus, 



N a(p) 2 ( N 1 N x k 

n 1 y 



and 



v{Pi) 2 
A 2 



(A21) 

(A22) 
(A23) 

(A24) 
(A25) 

(A26) 
(A27) 

(A28) 
(A29) 



Finally, the uncertainties of the linear ht parameters for fits to data with both correlated and random uncertainties 
(non-destructively readouts with read noise) are 



a 2 = <7 a (ran) 2 + a a (cor) 2 and 
of =a b (ran) 2 + a b (cor) 2 . 



(A30) 
(A31) 

The assumption that the uncertainties can be calculated separately for the correlated and random measurement 
uncertainties was tested via Monte Carlo simulations. Simulations for cases similar to that expected for the 70 /jm 
array are plotted in Figure IA16I As can be seen from these plots, equations IA30I and IA31I give very good estimates of 
the actual uncertainties. 
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